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LIST OF SYMBOLS AND ABBREVIATIONS 


vector quantities are denoted by arrows over the symbol 


ATP 

B 1 

b 


CRP 

Cp 

Gp 

dB 

dS 

f (x, t) = 0, f(y, t) = 0 


F(y ; x, t) 

H 

h 

n 

k(x, t) = 0, k(y, t) = 0 


K(y; x, t) 
M 


advanced technology propeller 

— ► 

i = 1,2; components of b in direction of principal 

curvatures 

| M t + | it r b = | b | 

b • n 

speed of sound in undisturbed fluid 

counter-rotation propeller 

power or torque coefficient of propeller 

thrust coefficient of propeller 

decibels 

element of blade surface area 

equation describing blade surface in reference 

frame fixed to the undisturbed medium 

f(y, t-r/c) = [f(y, t)] ret 

local mean curvature of blade surface 

IM n + IjCosq 

equation of surface defined such that f = 0 and 

k > 0 produces an open piece of the blade surface 

k(y, t-r/c) = [k(y, t)] fet 

local Mach vector based on c, M = IMI 

M • n 


I 


M 

M 

M 

M 


n 

r 


P 



NB 

OASPL 


P 

P 

P 

P B 0n> x) 


PCA 



r 

v 



SRP 


t 


time derivative of M 
M • n 
M- r 


projection of Mach vector on local normal plane 

to edges of blade surface, M = IM I 

projection of Mach vector on local tangent plane 

to blade surface, M^ = IM { I 

unit normal to blade surface 

number of blades 

overall sound pressure level 

acoustic pressure 

aerodynamic blade surface pressure 
time derivative of p 

blade surface pressure in blade fixed reference 
frame 

pitch change axis 
percentage chord coordinate 
x - y, r = lx - yl 
F/r 

/s -* 

r • v 

unit vector in the direction of the projection of r 
onto the local normal plane to the surface edges 
single-rotation propeller 
observer time 
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t 


1 



y 


Greek Symbols 

dy 

5(f) 

dZ 

0 

T\ 

^2 

A 

A 

A 0 


PO 

a b 


projection of r onto local plane tangent to blade 
surface 

forward flight velocity 

observer coordinate in ground fixed reference 
frame 

source coordinate in ground fixed reference frame 


element of length on edge of open surface 

Dirac delta function 

element of acoustic planform area 

— ► 

angle between n and F 

coordinate of blade surface point in blade fixed 
frame 

spanwise coordinate 
(1 + M n 2 - 2M n cos0) 1/2 
(A 2 + sin 2 0)l/ 2 

[Mp 2 cos 2 0 + (1 - Mp • FpSinY) 2 ]^ 2 

i = 1,2; components of in direction of principal 

curvatures 

inward geodesic normal of open surface defined 

by k>0 and f=0 

density of undisturbed fluid 

length parameter on f=0 along b 
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> t d 


°ir°22 


x 


Q. 

l 


K p K 2 

K M’ K r 


X 

x. 

0 ) 


two components of the tensor 
(tj t, - M ( M ( + tj M ( + M ( t, ) /A 2 
source time 

angle between r and edge of f=0 
n x o) 

principal curvatures of the surface f=0 
normal curvatures along M t , t^ , and b, 
respectively 
(cos0 - M )/A^ 

(cos0 + M n )/A 2 
angular velocity 
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INTRODUCTION 


In the late 1950’s and early 1960's the turbofan became the propul- 
sion device of choice for commercial aircraft fleets. The turbofan offered 
higher speed and a safer, more comfortable flight than its predecessor, the 
propeller. These conveniences, however, resulted in a substantial increase 
in fuel consumption. The uncertaintity of fuel supplies and the high cost of 
fuel in the 1970's revived interest in propeller driven flight. Researchers 
were confronted with the task of designing fuel efficient propulsion sys- 
tems that were capable of sustaining high speed and high altitude flight 
while maintaining acceptable ground and cabin noise levels. One outcome 
of this research was the development of the advanced technology propeller 
(ATP) or propfan. This new generation of propellers features highly loaded 
blades that are thin and swept back (Figure 1). It is estimated that the ATP 
will provide a 20 to 25% reduction in fuel consumption while maintaining 
the positive attributes of the turbofan [1]. A further reduction in fuel con- 
sumption can be achieved by eliminating the velocity swirl in the propfan 
wake. Counter-rotation propeller systems (Figure 2) accomplish this func- 
tion. Researchers estimate that, after acoustical treatment, the propfan can 
be made 8% more efficient by employing counter-rotation [2]. 

For the past several years considerable resources have been allo- 
cated by government and industry to produce a propfan design which 
minimizes noise while maintaining performance. The construction and 
acoustical testing of new propeller designs is an expensive and time 
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consuming process. Thus, a tractable, physically realistic analytical method 
for predicting propeller noise would provide the engineer with a cost effec- 
tive alternative to the above process. 

This work is concerned with the problem of predicting the noise gen- 
erated by advanced technology propellers operating at transonic/super- 
sonic helical tip speeds, in both the single-rotation and counter-rotation 
mode. The analytical method is based on two theoretical formulations 
developed by Farassat [3], [4]. One formulation is appropriate for subsonic 
sources and the other for transonic and supersonic sources. The computer 
program that implements these formulations, known as the Dunn-Farassat- 
Padula Advanced Technology Propeller (DFP-ATP) noise prediction pro- 
gram, was developed at the NASA Langley Research Center by the first 
author and F. Farassat and S. Padula of LaRC. Input to the program consists 
of blade kinematic and geometric quantities and either steady or periodi- 
cally unsteady blade surface pressure. Program output consists of pro- 
peller thrust and power coefficients, and the periodic acoustic pressure 
signatures and acoustic spectra for the thickness, loading and combined 
noise. 
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THEORETICAL FORMULATIONS 


Small perturbations in fluid quantities due to the motion of a body 
through an otherwise quiescent, ideal fluid are governed by the equations 
of linearized acoustics. Therefore, in the region exterior to the body the 

i 

acoustic pressure, p (x, t), is a solution of the homogeneous wave equation 
subject to the boundary condition that the fluid cannot penetrate the sur- 
face of the body. 

Using generalized function theory, the restricted domain of the above 
problem, i.e, the body's exterior, can be extended to include the entirety of 
space by the process of embedding [5]. Applying this procedure to the 
equations of acoustics leads to the Ffowcs Williams-Hawkins equation 
without the quadrupole term [6]: 

— 2 ' 

\r~2~ - V 2 p' = V [QlVflS(f)] (1) 

C at 

where 

Q = (-pff,M n ) , V 4 = (V,i-fj) 

and the bars over the differential operators denote generalized differenti- 
ation. The boundary condition and the pressure on the surface of the body 
appear as source terms in (1). The advantage of the embedding technique 
is that the Green's function for the wave equation in unbounded space is 


I 


7 


known, thus the solution of (1) can be obtained analytically provided the 
aerodynamic pressure on the surface of the body is known. 

A solution of (1), known as Formulation 1A, that is valid over source 
regions that move subsonically toward the observer, is given by [3]: 

p (x, t) = p T (x, t) + p L (x, t) 


where 


P T (x, t) = cp 0 


M (rM-?.+cM -cM ) 
n v i i r 7 


f=0 


2 3 

r (l-U/ 


dS 


Jret 


( 2 ) 


and 


P L < x - t) = c 


1_ 

/• 

PCOS0 

dS - 

r 

c 

J 


_r(l-M r ) Z 

ret ^ 

' L 


p(cos0-M ) 
r 2 (l-M r ) 2 


dS + 


Jret 


f=0 


f=0 


1 _ 

c 


r 


f=0 


pcos0(rM.7-+cM r *cM z ) 


2 3 

r (1-M/ 


dS 


Jret 


(3) 


The subscript ret indicates that the terms in the integrands are to be eval- 
uated at the emission (or retarded) time and the integrations are per- 
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I 1 

formed over the actual body surface. The functions p T and p L are referred 
to, respectively, as the thickness and loading noise components of the total 
pressure. This formulation is invalid for portions of the body that move 
toward the observer with transonic speed. A nearly equivalent formulation 
exists for supersonic sources, however, computational difficulties limit its 
usefulness. Particularly, instabilities in the numerical evaluation of (2) and 
(3) occur in regions where the Doppler factor, 1-M r , is "small", and it is 
computationally difficult to differentiate between "nearly" transonic and 
supersonic source regions. To circumvent these problems Farassat devel- 
oped equations (4) and (5) below [4]: 


2 Q m 

n N 


P T 


dX + 


^n ^F + ^F + ^F 


J 

F=0 

K>0 


Ar 


Ar 


dX - 


r M 2 Q c +M M 
n n av 


F=0 

K>0 


j 

F=0 

K=0 


V 


dy (4) 


where 


Pl = 


pQ 


■N 


fPQr 


dX + 


J 

F=0 

KM) 


Ar^ 


fpQi 


Ar 


dX + 


F=0 

KM) 


A 0 r 


dy 


F=0 

K=0 


(5) 


Q N =x[2>. 1 (-cos9-M n )+l] 


1 


Q R =- 

c 


U 2X 2 -^ 

A 


V 


M n + 2b K b + K l°ll + k 2°22 " 2Hh n + 


1 - 
-Q 

c 


M t - t l 

L X 2 


- 2Xb + 


f ^r+ 2XX. ' 
A 


J - 1 


Q c =2M - fXM - Q • b)+ k, S^+kJ^B 2 
r n|_c v n J 1 L 


9p 


+ TP d - b 


B 


c u 5o 


B 


q "=f (m - a . M t V K ™ M t 2 - 2HM n : 
F c^n t J m t n 


and 


Q c = XM + X, r 
av 1 \ 


This solution is referred to as Formulation 3. The surface integrals of 
equations (4) and (5) are evaluated over the acoustic planform (or 2-sur- 
face), and the line integrals are computed over the boundary of the acous- 
tic planform. The 2 -surface is generated by the intersection of the surface 
of the moving body and a sphere, centered at the observer position, x, 
whose radius collapses at the speed of sound. This formulation is valid in 
all regions regardless of the value of the radial Mach number. Due to the 
complex nature of the integrands, however, equations (2) and (3) can be 
evaluated much more efficiently in the subsonic regime. 
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per computes the acoustic pressure generated by advanced technology 
propellers. The contribution to the total noise attributable to subsonic por- 
tions of the propeller surface are computed via (2) and (3), while (4) and 
(5) are evaluated to determine the contribution from the transonic and su- 
personic regions of the propeller surface. 


COMPUTATIONAL STRATEGY 


In this section the computational strategies for the noise prediction 
process are discussed. The methods were designed to make the overall 
process as accurate as possible while minimizing the computational time 
and storage requirements. 

The calculations that are described here are based on two assump- 
tions. One, the observer position moves parallel to the propeller axis with 
the same speed as the aircraft, and two, for CRP predictions, both pro- 
pellers operate at the same RPM, they have the same number of blades, 
and their pitch change axes are initially aligned (Figure 3). The CRP re- 
strictions are imposed to facilitate the spectral analysis of the CRP acoustic 
pressure signal. If the user is interested in CRP predictions where these 
conditions cannot be satisfied, then the signals from the individual pro- 
pellers can be obtained by making two separate predictions, one for each 
propeller. The user can then add the two signals in an appropriate fashion 
and perform the desired Fourier analysis. In some instances the program 
can be easily modified to accomodate special circumstances. 

To begin the computations, user supplied blade geometry and loading 
data, which are given as functions of span and percentage chord, are mod- 
elled by high order, two dimensional, tensor product, least-squares splines. 
The software employed for this task is a modified version of the computer 
programs appearing in reference [7]. The blade surface is then separated 
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into two regions: the subsonic root region and the supersonic tip region. 

* 1 / 2 2\\I2 

The spanwise location = “fc* ' j » where c* = (l-£)c with 8 a 

small parameter supplied by the user, serves as the boundary separating 
the two regions. This separation point avoids the computational difficulties 
that were discussed in the previous section. Both regions are then subdi- 
vided into equally spaced spanwise and chordwise intervals according to 
user specifications (Figure 4a). For each resultant panel, the following cal- 
culations are performed: 

STEP (1 ) A user-specified number of Gauss-Legendre (GL) nodes are dis- 
tributed over the panel (Figure 4b). 

STEPT2) Blade geometric quantities and the steady portion of the blade 
surface pressure at each GL node are computed and stored. 

STEP(3^ For each observer time, the emission time(s) and are 

calculated at each GL node. If 1-M f < 8 on any GL node at any 
emission time, then STEP(3A) is performed, otherwise STEP(3B) is 
performed. 

STEP(3 A! The panel is subdivided further, according to user specifications, 
by evenly spaced spanwise and chordwise partitions (Figure 4c) and 
equations (4) and (5) are used to compute the noise for that 
observer time. The surface integrals are evaluated by the "collapsing 
sphere" method [8]. The T -curves generated by the method are 
approximated by piecewise linear functions and the integrals are 
approximated by the mid-point rule. The line integrals appearing in 
(4) and (5) are calculated by the trapezoid rule. 

STEPf3B) Equations (2) and (3) are used to compute the noise in this case. 
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The integrals are calculated by Gaussian quadrature based on the 
nodes computed in STEP(l). The software used to calculate the 
Gaussian quadrature nodes and weights was taken from [9]. 

STEPT41 For each observer time, the noise contribution from each panel is 
summed to form a single blade pressure signature. This is done for 
the thickness, loading and total noise components. For CRP 
predictions the above 4 steps are performed twice, once for each 
propeller. The single blade signatures of the two propellers are then 
summed according to the superposition principle. 

STEP(51 The single blade signatures are shifted and summed to obtain an 
NB blade periodic acoustic pressure signal, again according to the 
superposition principle. 

STEP(61 The Fourier series of the NB blade signal is computed via Simp- 
son's rule. 

STEP(71 Sound pressure levels are computed and written to output as are 
the NB blade thickness, loading and combined acoustic pressures. 

The above steps are depicted schematically in Figure 5. As was men- 
tioned earlier, the integrands of equations (2) and (3) must be evaluated at 
the emission time. Also, construction of the L-surface of equations (4) and 
(5) is based on a knowledge of the emission times. Thus, it is important to 
be able to calculate the emission times efficiently. The algorithm employed 
for this purpose is responsible for approximately one-half of the program's 
total computational time and is briefly described in appendix A. 
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PROGRAM DESCRIPTION 


In this section the user input specifications, program output and de- 
tails of the program usage are described. 


Operating Instructions 

The user input file is associated with the FORTRAN logical unit num- 
ber 5. All input parameters reside on this file. User input consists of three 
groups of NAMELIST variables and blade coordinates and loading data. The 
structure and contents of the input file will be discussed in the ensuing 
subsections. The standard program output file is associated with the 
FORTRAN logical unit number 6. Additional output files containing inter- 
mediate results are associated with FORTRAN logical unit numbers 7 and 8. 
The user must open these files in the calling program. 

To access the noise prediction program, the user's calling program 
must have the FORTRAN statement 

CALL DFPATP( X , NXDIM ) 

where X is a one-dimensional dynamic storage array whose dimension, 
NXDIM, must be at least 

max{ 3*NOBS + 6*NT + 26*N1S*N2S + 3*max(NlS,N2S) + 7*(IGAUSS) 2 , 

6*NEPTS *(NQPTS+ 1 ) + 2*NQPTS } . 
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The standard program output consists of a banner page, a one line 
identifier, NAMELIST parameters, computational grid information, power 
(or torque) and thrust coefficients, and thickness, loading and combined 
sound pressure levels and periodic acoustic pressure signatures. In addi- 
tion, unit 7 contains the single blade pressure signal for SRP predictions. 
For CRP predictions, unit 7 contains the single blade signal for the forward 
propeller and unit 8 contains the single blade signal for the aft propeller. 

NAMELIST PHYSCAL Parameters 


The variables appearing in this NAMELIST describe the physical 
characteristics of the undisturbed medium and the kinematic properties of 
the propeller. 


Variable 

BTIP 

B INNER 

VF 

RPM 

CO 


Description 

Real. Distance from propeller hub to propeller tip in 
meters. 

Real. Distance from propeller hub to propeller root in 
meters. 

Real. Propeller forward flight velocity in meters per 
second. 

Real. Propeller angular velocity in revolutions per minute. 
Real. Speed of sound of undisturbed fluid in meters per 
second. 



RHOO 

NB 

NOBS 

STEADY 

NCRP 

PROPDIS 


Real. Density of undisturbed fluid in kilograms per cubic 
meter. 

Integer. Number of blades. 

Integer. Number of observers. 

Logical. STEADY = .TRUE. => steady loading. 

STEADY = .FALSE. => unsteady loading. 

Integer. NCRP = 0 => SRP prediction. 

NCRP = 1 => front propeller of CRP prediction. 

NCRP = 2 => aft propeller of CRP prediction. 

Real. Distance between CRP propellers in meters 
(Figure 3). 


NAMELIST GRID Parameters 

These parameters determine the computational grid for the acoustic 
calculations. 

Variable Description 

NT Integer. Number of observer time points for single blade 

pressure signature. NT > 2. Default value: 125*NB. 

NS Integer. Number of spectral levels to be output. Default 

value: 30. 
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N1SUB 

N1SUP 

N2SUB 

N2SUP 

N1S 

N2S 

EPSILON 

IGAUSS 


Integer. Number of equally spaced spanwise intervals for 
the subsonic portion of the blade (Figure 4a). Default 
value: 2. 

Integer. Number of equally spaced spanwise intervals 
for the supersonic portion of the blade (Figure 4a). 

Default value: 2. 

Integer. Number of equally spaced percentage chord 
intervals for the subsonic portion of the blade (Figure 4a). 
Default value: 4. 

Integer. Number of equally spaced percentage chord 
intervals for the supersonic portion of the blade 
(Figure 4a). Default value: 4. 

Integer. Number of equally spaced spanwise points for 
Formulation 3 panels (Figure 4c). Default value: 10. 
Integer. Number of equally spaced percentage chord 
points for Formulation 3 panels (Figure 4c). Default 
value: 10. 

Real. If EPSILON > 1-M r for some point on a panel then 
Formulation 3 is used. Default value: 0.05. 

Integer. Order of Gaussian quadrature for computing 
Formulation 1A integrals (Figure 4b). Default value: 7. 


NAMELIST OBSERV Parameters 


This NAMELIST contains the observer(s) coordinates. 
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OBS(I, J) Real. I-th cartesian coordinate, in meters, of the J-th 

observer with respect to the aircraft-fixed reference 
frame (see below). I = 1,2,3 and J = l,...,NOBS. 

Blade Loading and Geometry Input Specifications 

Details of the blade geometry input begin with a discussion of the co- 
ordinate systems used in the acoustic calculations. They are the ground 
fixed x-frame, the aircraft fixed X -frame and the blade fixed rj-frame 
(Figure 12). Initially, the propeller hub center is at the origin of the 
x-frame. The x^ axis and the propeller axis coincide, with positive x^ in the 
direction of flight. The propeller disk lies in the x j _x 2 P^ ane initially and 
Xj and *2 are defined such that the x-frame is right-handed. The positive 
sense of rotation is defined to be from the x^-axis to the X 2 -axis. The 
X-frame is identical to the x-frame initially and moves in the forward 

flight direction with speed Vp. The observer coordinates are given in the 

— > _ 

X-frame. The origin of the blade fixed Ti-frame coincides with the origin of 
the X-frame at all times. The ri^-axis lies along the propeller axis, the 
r^-axis lies along the pitch change axis, and is chosen so that the 
rj-frame is right-handed. 

The blade geometry input is specified with respect to the blade fixed 
frame. Airfoil sections are generated by cutting the PCA with planes per- 
pendicular to the PCA. For each section the twist, a, in degrees, the local 
chord length, CHORD, in meters, and the rjj and coordinates of the lead- 
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ing edge, rj^-LE and rj^-LE, in meters, are specified (Figure 13a and 13b). 
The cutting planes are given by 11 ^= constant (meters). The first cutting 
plane is = BINNER and the last is = BTIP. The coordinates of the air- 
foil section are nondimensionalized by the local chord length. The upper 
and lower coordinates of the airfoil, y^(Q) and y^(Q), are given as functions 
of percentage chord, Q, with Q = 0 corresponding to the leading edge and 
Q = 1 corresponding to the trailing edge (Figure 13c). 

Blade surface pressure can be either steady or periodically unsteady. 
The user must supply the steady portion of the surface pressure at each 
(^ 2 * Q) point that is input. For each spanwise cut the upper and lower 
steady surface pressure, py(Q) and 

2 2 

centage chord. These quantities are nondimensionalized by PqCO BTIP . 

If the surface pressure is unsteady, then it must be of the form 
P(^2’ Q)[l + SCOl- The function g is assumed to be periodic with period 2k. 
The user specifies g in tabular form, (t., g(t.)), with the first time point 
equal to 0 (zero) and the last equal to 2k. 

The geometry and loading portion of the input file has the following 
format: 


p^(Q), are input as functions of per- 


(1) Number of = constant cutting planes (NEPTS), number of Q points 
per airfoil section (NQPTS), number of time points for unsteady 
loading (NTPTS). - FORMAT(1X,3I10) - 

The chordwise Q points should be biased toward the leading edge 
and the spanwise points should be biased toward the tip. 
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(2) Percentage chord grid. — F0RMAT(1X,E15.7) 

Q x = 0.0 

q nqpts = 10 

(3) r| 2 , CHORD, a, rij-LE, ti 3 -LE - F0RMAT(1X,5E15.7) 
y L (Q iX y^Qj), P^Qp, P^Qp -- FORMAT(ix,4Ei5.7) 

y L^NQPTS^’ y U^NQPTS^’ P L^NQPTS^ P U^NQPTS^ 

(4) Item (3) is repeated for each of the NEPTS spanwise points. 

(5) If the surface pressure is unsteady, then g(t) is input here. 

t r g(tj) - FORMAT( 1X,2E1 5.7) 


^TPTS’ g(t NTPTS ) 
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( 6 ) 


If the noise prediction is for a CRP configuration, then items (1) to (5) 
are repeated for the aft propeller. The location of the aft propeller 
data within the input file is described in the Input File Structure 
subsection. 

The geometric and loading input parameters are summarized below. 


Variable 

NEPTS 

NQPTS 

NTPTS 

Q 

^2 

CHORD 

a 

t \ 1 - LE 

„ 3 -le 

y L (Q) 


Description 

Integer. Number of spanwise data points. 

Recommend: 20-30. 

Integer. Number of percentage chord data points 
Recommend: 30-50. 

Integer. Number of time points for unsteady loading data. 
Recommend: 20-50. 

Real. Percentage chord data points. Nondimensionalized 
by CHORD. 0 < Q < 1. 

Real. Spanwise data points in meters. 

BINNER < T] 2 < BTIP. 

Real. Local chord length in meters. 

Real. Local twist angle in degrees. 

Real, rij coordinate of leading edge in meters. 

Real, ri^ coordinate of leading edge in meters. 

Real. Coordinate of lower airfoil surface. 
Nondimensionalized by CHORD. 
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y Tj(Q) 

P L (Q) 

Pu(Q) 

t 

g(t) 


Real. Coordinate of upper airfoil surface. 
Nondimensionalized by CHORD. 

Real. Steady pressure component on lower airfoil surface. 

2 2 

Nondimensionalized by PqBTIP to . 

Real. Steady pressure component on upper airfoil surface. 

2 2 

Nondimensionalized by p^BTIP go . 

Real. Time coordinate for unsteady portion of blade 
pressure. Nondimensionalized by 27 c/cd. 0 < t < 2 jc . 

Real. Unsteady portion of blade pressure. Nondimensional. 


Input File Structure 

The structure of the user input file is described here. Two sample input 
files are presented in the next section. 

(A) SRP predictions 

(1) An 80 character identifier with A80 format 

(2) NAMELIST PHYSCAL - Note that NCRP must be equal to 0. 

(3) NAMELIST GRID 

(4) Blade coordinates and steady portion of the blade loading data 
(4') If STEADY = .FALSE., then the unsteady portion of the blade 

loading data appears here. 

(5) NAMELIST OBSERV 
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(B) CRP predictions 


(1) An 80 character identifier with A80 format 

(2) PHYSCAL for forward propeller. — NCRP must be equal to 1. 

(3) GRID for forward propeller 

(4) Forward propeller blade coordinates and steady portion of the 

blade loading data 

(4’) If STEADY = .FALSE., then the unsteady portion of the forward 
propeller blade loading data appears here. 

(5) OBSERV — The observer(s) coordinates are based on the hub of 

the forward propeller being at the origin of the aircraft fixed 
frame (Figure 3). 

(6) PHYSCAL for aft propeller -- NB and NOBS must be the same for 

both propellers, NCRP must be equal to 2, and RPM for the aft 
propeller must be the negative of the forward propeller RPM. 

(7) GRID for aft propeller — NT must be the same for both 

propellers. 

(8) Aft propeller blade coordinates and steady portion of the blade 

loading data 

(8') If STEADY = .FALSE., then the unsteady portion of the aft 
propeller blade loading data appears here. 

(9) OBSERV — These data should be identical to the data for the 

forward propeller. 
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EXAMPLES 


(A) SRP prediction 

To illustrate the SRP prediction capability, the noise generated by a 
model SR3 propeller was computed. The steady loading used in the calcu- 
lations was obtained from Denton's code [10]. The computed power and 
thrust coefficients were Cp = 1.068 and Gp = 0.250. The overall sound 
pressure level was OASPL = 142.15 dB. Graphs of the sound pressure level 
harmonics and periodic acoustic pressure signals for the thickness, loading 
and total noise components can be seen in Figures 6, 7, and 8. The input 
stream for this example is shown below. 

SRP PREDICTION --SR3 PROPELLER 

$PHYSCAL 

BTIP = 0.312, 

BINNER = 0.0936, 

VF = 241.2, 

RPM = 7878.4, 

CO = 310.5, 

RHOO = 0.457, 

NB =4, 

NOBS = 1, 

STEADY = . TRUE., 

NCRP = 0, 

PROPDIS = 0.0, 

SEND 
$GRID 
NT = 500, 
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NS = 30, 

N1SUB = 2, 

N1SUP = 3, 

N2SUB = 4, 

N2SUP = 4, 

N1S = 10, 

N2S = 10, 

EPSILON = 0.05, 

IGAUSS = 7, 

$END 

• —blade geometry 

• and loading 

• data 

SOBSERV 
OBS(l,l) = 0.810, 
OBS(2,l) = 0.0, 

OBS(3,l) = -0.20, 

$END 


(B) CRP prediction 

These results were generated by a model advanced counter rotation 
propeller system. The blade loading, which includes the aerodynamic in- 
teraction between the forward and aft propeller, was obtained from 
Adamczyk’s Euler code [11]. For the forward propeller Cp = 1.225 and 
Crj, = 0.332, and for the aft propeller Cp = 1.296 and Gp = 0.375. The OASPL 
was 146.33 dB. Graphs of the noise components appear in Figures 9, 10, 
and 11. The input file is presented below. 
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CRP PREDICTION 


$PHYSCAL 

SPHYSCAL 

BTIP = 0.3064, 

BTIP = 0.2937, 

BINNER = 0.137, 

BINNER = 0.1322, 

VF = 242.9, 

VF = 242.9, 

RPM = 8283, 

RPM = -8283, 

CO = 338.3, 

CO = 338.3, 

RHOO = 1.039, 

RHOO = 1.039, 

NB = 8, 

NB = 8, 

NOBS = 1, 

NOBS = 1, 

STEADY = . TRUE., 

STEADY = .TRUE., 

NCRP = 1, 

NCRP = 2, 

PROPDIS = 0.1037, 

PROPDIS = 0.1037, 

$END 

SEND 

$GRID 

SGRID 

NT = 1000, 

NT = 1000, 

NS = 30, 

NS = 30, 

N1SUB = 2, 

N1SUB = 2, 

N1SUP = 2, 

N1SUP = 2, 

N2SUB = 4, 

N2SUB = 4, 

N2SUP = 4, 

N2SUP = 4, 

N1S = 10, 

N1S = 10, 

N2S = 10, 

N2S = 10, 

EPSILON = 0.05, 

EPSILON = 0.05, 

IGAUSS = 7, 

IGAUSS = 7, 

$END 

SEND 

• forward blade 

• ---aft blade 

• geometry and 

• geometry a 

• loading data 

• loading dai 

SOBSERV 

SOBSERV 

OBS(l,l) = 0.800, 

OBS(l,l) = 0.800, 

OBS(2,l) = 0.0, 

OBS(2,l) = 0.0, 

OBS(3,l) = -0.05, 

OBS(3,l) = -0.05, 

SEND 

SEND 

continued on next column 



27 


DISCUSSION 


The effects of varying the computational grid parameters N1SUB, 
N1SUP, N2SUB, N2SUP, IGAUSS, N1S, N2S and EPSILON have been studied 
and the results reported in reference [12]. Though the study was not ex- 
haustive, it was determined that acceptable accuracy levels could be at- 
tained without the consumption of large amounts of computer resources by 
using the default values listed in the Program Description section. Fur- 
thermore, the authors speculate that a comprehensive grid study will not 
improve the efficiency of the program to an appreciable degree. To some 

extent, the choice of these parameters is problem dependent. For example, 

* 

in situations for which the "sonic line", , is closer to the root of the 

propeller than the tip, the user would be justified in adding one more 
spanwise interval to the supersonic region of the blade, i.e., set N1SUP 
equal to 3, and for strictly subsonic predictions adequate results may be 
obtained by reducing IGAUSS, which would result in a substantial decrease 
in computational time. However, for most propfans operating at their 
cruise conditions, the default parameters will yield satisfactory noise 
predictions. 

Comparisons between the theoretical formulations described here 
and experimental data, both wind tunnel and actual flight data, have been 
conducted [12]. Measured acoustic data, however, usually contain the ef- 
fects of nearby hard surface reflection, wave propagation in the fuselage 
boundary layer or other physical phenomena as well as additional sources 
of noise such as engines. To adequately model the noise produced in situa- 
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tions where these effects are important, the noise prediction program 
should be used in conjunction with other analytical or empirical techniques 
that sufficiently describe the above physical aspects. 

In summary, an efficient computational method for predicting the 
noise produced by high speed, advanced technology propellers is pre- 
sented. This new generation of propellers features thin and highly swept 
blades. Therefore, the noise generated by the propfan can be adequately 
modelled by linear acoustics. Two analytical solutions of the equations of 
linear acoustics have been developed by Farassat. One solution. Formula- 
tion 3, describes the noise produced by transonic and supersonic propeller 
motion and the other. Formulation 1A, describes the noise generated by 
subsonic propeller motion. The program decides automatically which of the 
two solutions to calculate. The calculations were designed to make the 
evaluation of the integrals appearing in equations (2) to (5) as accurate as 
possible while minimizing CPU resources and computer memory require- 
ments. In addition, efforts were made to keep user interaction simple and 
program output descriptive. 
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APPENDIX A Emission Time Calculation 


For a given observer time, t, consider a particular source point and 
observer point pair. If sound is emitted from the source at time x and ar- 
rives at the observer at time t then t is called the emission time (or re- 
tarded time). Let y(x) and x(t) be the cartesian coordinates of the source 
point and the observer point, respectively, then, since the motion of the 
source and observer are known, the emission time can be determined from 
the equation 

c( t - x ) = I x(t) - y(x) I . (6) 

For propeller motion, equation (6) can be written as 

A<|> + B<|) + C + cos( <j> + D) = 0 , (7 ) 

where <J> = co( x - t ) and the coefficients A, B, C, and D are known functions 
of the source and observer positions and the parameters of motion. Equa- 
tion (6) or (7) is referred to as the "retarded time equation". Note that for 
propeller motion the emission time appears as an implicit variable in the 
retarded time equation. Equation (7) cannot be inverted analytically, 
therefore the emission times must be determined through a numerical 
technique. 

It was previously mentioned that for every observer time compris- 
ing the acoustic pressure period, the emission time(s) must be calculated 
for each source point on the propeller blade computational grid. For typical 
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noise predictions, the retarded time equation must be solved on the order 
6 7 

of 10 to 10 7 times accounting for nearly 50 percent of the program's total 
computational time. The algorithm outlined below was designed to exhaus- 
tively search for the solutions of the retarded time equation in an efficient 
manner. 

2 

Let h(4>) = A<J> + B<j) + C + cos( <|> + D) . Then, solutions of h(<|>) = 0 are 

<l> 

sought such that — < 0 . Due to causality, only solutions of this type are 
physically plausible. The solutions of h = 0 are found by a hybrid Bisec- 
tion/Newton method. Since multiple solutions are possible, it is necessary 
to obtain "close" initial estimates for each root so that Newton’s method can 
effectively locate all of the solutions. This is accomplished by finding inter- 
vals in which the curvature of h does not change sign. Within such inter- 
vals, h exhibits parabolic behavior. The function is then sampled at the 
endpoints of an interval, and several tests are conducted to determine if 
there are zero, one or two roots in the interval. If roots are detected in the 
interval, then h is approximated by a parabola and the root(s) of the 

parabola is (are) used as the initial estimate(s) for Newton's method. New- 

dh 

ton’s method can be unstable if “ is "small" in some neighborhood of the 
root, thus an alternative method is required for these situations. The Bi- 
section method is used to further refine the estimate before using New- 

d h 

ton's method when ~ is "small". This procedure is repeated for each inter- 

d(p 

val in which the curvature does not change sign. 
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Figure 2: Counter Rotation Propeller System 
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Figure 3: Counter Rotation Configuration - Aircraft Fixed Frame 
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Figure 4: DFP-ATP Computational Grid 
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Figure 5: DFP-ATP Flow Chart 
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Figure 8: Total Noise for Example 1 - 

Model SR3 Propeller 




Figure 9: Thickness Noise for Example 2 - 

Model Counter Rotation Propeller System 
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